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BORIS KRUGLIKOV 

Abstract. In this paper we investigate compatible overdeter- 
mined systems of PDEs on the plane with one common charac- 
teristic. Lie's theorem states that its integration is equivalent to a 
system of ODEs, and we relate this to the geometry of rank 2 dis- 
tributions. We find a criterion for integration in quadratures and 
in closed form, as in the method of Darboux, and discuss nonlinear 
Laplace transformations and symmetric PDE models. 



Introduction 

Consider a system £ of partial differential equations of (maximal) 
order k (geometrically a submanifold in fc-jets). We will study its inte- 
gration micro-locally near a regular point (in a neighborhood U d £). 

We will mainly restrict to the case of two independent and one de- 
pendent variables (non-scalar systems can be treated similarly or re- 
written via Drah's trick [StoJ; the case of more independent variables 
makes no fundamental difference either). 

0.1. Formulation of the problem. We assume the system £^ is overde- 
termined and is compatible (formally integrable). 

Whenever the above assumptions on £ are fulfilled we can evaluate 
formal dimension and rank of the system (like in Cartan test IC3I ), 
namely determine on how many functions of how many variables a 
general solution formally depends. 

For a (determined or overdetermined) system on the plane two situ- 
ations are possible. Either the system has empty characteristic variety 
or the latter is an effective divisor on CP^ (Weil divisor in algebraic 
geometry; here a finite set of points with multiplicities). 

In the first case the system £ is of finite type, i.e. its local solution 
space is finite dimensional. Then integration of £ can be reduced to a 
system of ODEs; we call such systems Frobenius. 



Key words and phrases. Lie's class 1, Darboux integrability, system of PDEs, 
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In the second case additional assumptions must be imposed to guar- 
antee existence of solutions. Counting of solutions can still be carried 
on the formal level and the general stratum of the solution space of £ 
is parametrized by uj functions of 1 varia blfl where oj is the degree of 
the divisor, i.e. the number of points counted with multiphcities. 

Note that u can be described in a different way: Since the charac- 
teristic variety is discrete, the symbol gu of the system £ stabilizes and 
its dimension is dim^f^ = u for big k. 

In this paper we restrict to the case w = 1, which is the next simplest 
after the finite type case w = 0. 

0.2. Main results. For uj = 1 Sophus Lie obtained in 1893 a theorem, 
which states that this case can be reduced to ODEs as well: 

Theorem A. A compatible regular overdetermined system £ of class 
a; = 1 can be locally integrated via ODEs. 

The proof in [L] is rather sketchy. The result was later obtained in 
|De] without reference to Lie. We will demonstrate the claim via geom- 
etry of rank 2 distributions and relate it to other important result^. 

Three remarks are of order. First is that we can change here formal 
integrability to local integrability. Generally it is wrong due to Levi 
and similar examples. The essential feature is regularity of the char- 
acteristic variety and peculiarity of a; = 1 (i.e. for a; > 1 the passage 
from formal integrability to local integrability is not generally correct). 

Second, and more important, is that the reduction procedure can be 
made explicit and this allows not only to claim the reduction theoreti- 
cally, but also to develop a practical algorithm for integrability. 

And third is that facing only the problem of reduction, we can allow 
arbitrary number of independent and dependent variables, leaving only 
the requirements of regularity and a; = 1, i.e. precisely one common 
characteristic (counted with multiplicity) for the equations of £. 

In this paper we mostly concentrate on the problem of effective inte- 
gration of PDE systems. So we are especially interested in quadrature 
and integration in closed form (an ideal case of Darboux integrability 
- the definitions are below in the text, see also |Dat [F| lAFVj ). 

Theorem B. A system £ of class u = 1 is integrable in closed form 
iff it is linearizable by an internal transformation. 

^Sophus Lie called the number w class of the system; for EUie Cartan this is the 
character Si (provided the Cartan number is 1: = for i > 1). 

^Let us mention that paper ^ discusses another reduction to ODEs for the 
involutive PDE systems of the 2nd order. This family meets ours by class ui = 1 
systems of type 2E2 in terminology of W2.l[ which were studied by Cartan in |Ci| . 
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A system £ of class 00 = 1 is integrable in quadratures iff it has a 
transitive solvable Lie algebra of internal symmetries^ 

Of course, one is interested in algorithmic integration, so that an 
effective hnearization is important. Then a sequence of generahzed 
Laplace transformations (these are the external transformations intro- 
duced for class w = 1 in p^]) finishes the job. 

As we shall explain, determining both linearization and quadrature 
is related to investigation of the rank 2 distributions internally related 
to the system £. Linearizable systems correspond to Goursat distribu- 
tions, i.e. canonical Cartan distributions on the jet spaces for ODEs (in 
general non-linear situation the growth vector is un-restricted) . Rank 2 
distributions for the systems integrable in quadratures have the struc- 
ture of integrable extensions, which can be decoded starting from its 
Tanaka algebra. 

Thus we can model types of reduction, based on the normal forms of 
rank 2 distributions. In particular, the simplest among exactly solvable 
class uj = 1 compatible non-linearizable PDE systems will be those that 
can be reduced to Hilbert-Cartan equation (it's symmetry algebra will 
contain the exceptional Lie group 6*2)- More complicated examples 
will be presented at the end of the paper. 

0.3. Structure of the paper. We will exploit the geometric theory 
of PDE, jet-geometry and the basics of Spencer formal theory. We will 
also use the geometry of vector distributions. The reader is invited to 
consult [H [KLV] for details. 

Notations are different from source to source, and we adapt those of 
[KL] . Since this paper is a continuation of p^], an acquaintance with 
the latter will be useful (but not mandatory). 

The paper is organized as follows. 

In Section [1] we recast the class uj = 1 systems into the language 
of the geometry of differential equations and provide a new modern 
proof of Theorem \^ Reduction to rank 2 distributions is the crucial 
ingredienlQ. We then discuss an algorithmic method to integrate such 
systems and prove the 1st half of Theorem [Bl 

In Section [2] we discuss another more general method of integration 
of PDEs via integrable extensions (coverings), and relate this to the 

■^Internal symmetries are transformations of the equation E considered as a man- 
ifold preserving the induced Cartan distribution. They are more general than the 
classical Lie symmetries, but can differ from the higher (Lie-Backlund) symmetries. 

^There is no difference in arguments and in Section [T] we take dimension of the 
base K" {x) arbitrary. In further sections we restrict to n = 2 independent variables 
for simplicity of exposition. 
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generalized symmetries. Notice that integrable extensions for rank 2 
distributions were classified in |AKj . so their description in the sym- 
metric cases reduce to purely algebraic questions. 

In Section [3] we formulate the main invariants of compatible systems 
£ of class w = 1, and we discuss transformations of such systems in 
linear and non-linear cases. We investigate linear system from the 
viewpoint of internal geometry (complimentary to the external point 
of view in [K]), obtain the linearization criterion and finish the proof of 
Theorem [Bl Depending on the type of the system and its reduced rank 
2 distribution we can describe the structure of the general solution and 
a method of its integration. 

Section m is devoted to various examples of compatible PDE systems 
of class uj = 1. We will perform integration via the method of integrable 
extensions, generalized nonlinear Laplace transformations and discuss 
their relation to Darboux integrability. Some of the most symmetric 
examples are coverings of the overdetermined involutive system of order 
2 on the plane investigated by E.Cartan. 

Acknowledgment. It is a pleasure to thank Nail Ibragimov for his 
translation of Sophus Lie paper [L] in |LGj , which was a starting point 
for this paper. I am grateful to Valentin Lychagin for many discussions 
on the initial stage of the project. The paper is strongly influenced by 
a collaboration with Ian Anderson, to whom I am thankful. 

Hospitality of MLI (Stockholm) in 2007, MFO (Oberwolfach) in 
2008, Banach center (Warsaw) and Utah State University in 2009, and 
IHES (Paris) in 2010 is highly acknowledged. 

1. Around Sophus Lie theorems 

In this section we give a modern proof of Theorem |X1 Sophus Lie's 
original approach is indirect and hard to implement. We present a 
geometric method, which is the base of our approach to integration of 
class a; = 1 systems. Furthermore we will elaborate this theorem to 
get the constructive Theorem [Bl 

1.1. The geometric setup. Consider a compatible overdetermined 
system of PDEs as a submanifold in the space of jets £ C J^{W,N), 
where = is the space of independent variables x = (x*), = 
is the space of dependent variables u = {u^) (assuming £ to be of pure 
order k is not crucial). Let vr^ : J'^ — )■ W , T^k,k-i '■ J^~^ denote 

the natural projections. 

We prolong the system to the level it becomes involutive (see the 
discussion about relation of this with compatibility in |KL] and [K]). 
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The assumption on the class a; = 1 yields dimgk = 1 for the symbol 
space starting from this level (still denoted by k). 

This jet-space is equipped with the canonical Cartan distribution C = 
Ann(^^ : \a\ < k) C. TJ^ (where 6^. = du{. — Y^u-'^+^Ax'^ in canonical 
coordinates). The induced Cartan distribution on the equation Cs = 
C n TS has rank n + 1. Indeed, it is generated by one vertical vector - 
a generator of the symbol gk = Kei^diik^k-i ■ TS — )■ TJ^^^), which is 
defined up to scaling, and n total derivatives 'Dx^ = d^-i + 'Y'^i+u^u^ 
restricted to the equation, which are defined mod Qk- Denoting the 
horizontal space by H (not canonical) we get 

Ce = H® gk. 

Lemma 1. Let p G PT*W be the (unique) characteristic covector. 
There is a unique (n — 1)- dimensional subdistribution U G C£ such that 
duklll) = p"*" C TW and U consists of Cauchy characteristics of Cg. 

Proof. Let H be some choice of horizontal space, Tl d H the lift of p-^ 
and rj a vertical vector field (section of gk). 

Let us use the standard identification 7r^j._i(*) — S^T*W ® TN 
[KEj. Then the condit ion (jj = 1 translates to gu = {p^ ® v) for some 
V e TN. The Lie bracket induces a pointwise bracket H^g^ — > i^k-i = 
(/-I 0v) C S^-^T*W ® TN. 

With identification H ~ TW this latter is the restriction of the 
natural pairing TW O S^T*W O TiV ^ S''-^T*W O TN. It follows 
that [^,r]] = modCe and = modC^ V G U. 

It remains to choose an additional vector field in if \ 11 and consider 
the induced bracket r : (C) (g) 11 — )■ Uk-i- It can be non-zero since the 
pairing (C,P^ ® 9) kp'^~^p{() ^6^0. 

Let us change H = (C) © 11 by modifying 11 as the graph of the map 
— r(C, •) e n*®t'fc-i — n*(8)(?fc (identification via -Z^J^;,). Then the new 
space n is still involutive with respect to the induced bracket and it 
commutes with both ( and r] mod C^. This means that the sections of 
n are Cauchy characteristics. 

Uniqueness of 11 follows from the fact that the above (bracket) pair- 
ing {() ® gk ^ ^k-i is non-zero. □ 

Remark 1. For n = 2 characteristic vectors are dual to characteristic 
covectors. It is not however true that the former can be lifted to Cauchy 
characteristics of C^. This is peculiarity of the case cu = 1. 

1.2. Reduction to rank 2 distributions. Due to Lemma [1] internal 
geometry of the distribution Cs is equivalent to that of the rank 2 
distribution C^/II. This implies Sophus Lie theorem: 
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Proof of Theorem [Al Consider the pair (£,Cf). Solutions of 
the system are n-dimensional integral submanifolds of the distribution, 
whose projection to the base are submersive. 

It follows from the proof of Lemma [1] that an n-dimensional subspace 
of Cf is involutive with respect to the (bracket) pairing A^C^ — )■ Uk-i 
iff it contains 11. In other words, a solution must be tangent to 11. 

It is the standard fact, that the sub-distribution 11 generated by 
Cauchy characteristics is integrable and shifts along it are symmetries 
for Cs- Taking the (local) quotient we arrive to the manifold M = S/U 
(quotient by the leaves) equipped with a rank 2 distribution A = Cs/U 
without characteristics. 

Such a distribution has integral curves, which are found by solv- 
ing underdetermined ODEs. The space of integral curves is locally 
parametrized by 1 arbitrary function of 1 variable (determinization of 
the ODE, given e.g. by the choice of a curve in the image of any sub- 
mersion M — > with fibers transversal to A). The inverse of the 
quotient map £ ^ M sends any of them to an ra-dimensional integral 
surface, i.e. the solutions of £. □ 

Note that this proof, as well as the arguments from the previous 
subsection, uses integration of ODE systems twice: first to solve the 
Frobenius system, corresponding to Cauchy characteristics 11, and then 
to find the integral curves of A = Cf/II. 

The latter integration can be split in turn into integration of the 
bracket-closure of the distribution Aoo = A + [A, A] + . . . , which is 
Frobenius and then integrating A in the leaves. 

In the first case the order of the system is dim£^ — (^ — 1) = codim 11. 
In the second it is split into an ODE of order equal to the number of first 
integrals for A in M (= codim Ago) and an ODE of order dim A^o — 2. 

1.3. Constructive integration methods. A theorem of Sophus Lie 
states that ODEs with a transitive solvable Lie algebra of symmetries 
are integrable in quadratures. This is equivalent to the claim that if a 
holonomic distributiorH A on a manifold M admits a solvable symmetry 
Lie group of complimentary dimension with orbits transversal to it, 
then the integral leaves of A can be expressed in quadratures [KLR] . 

We extend this theorem to non-holonomic distributions A. We as- 
sume at first the distributions are completely non-holonomic, i.e. the 
bracket closure Aqo equals TM and so A has no first integrals. Generic 
such distributions have no integral surfaces, and integral curves (which 
always exist) are the maximal integral manifolds. 



^This means it satisfies the Frobenius condition [r(A),r(A)] C r(A). 
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Theorem 1. Let A be a completely non-holonomic distribution of rank 
r on a manifold M"- . Suppose a solvable Lie group G of dimension n — r 
acts by symmetries with orbits everywhere transversal to A. Then local 
integral curves of A can be found by quadratures. 

Proof. Denote n : M ^ = M/G the local quotient by the orbits, 
(the space of G- invariants). Notice that vr* maps A to TL. 

Choose a curve 7 C L and restrict the distribution A to 7r^^(7). 
This is a line field and G acts transitively by symmetries on 71^^(7). 
By the classical Sophus Lie theorem the integral curves of this line field 
can be found by quadratures. Thus these restricted integral curves are 
parametrized by n — r integration constants in 7r~^(7), while the curves 
7 C L are parametrized by r — 1 function of 1 variable. The integral 
curves of A in M are given through these by quadratures. □ 

In particular, for our case r = 2 we get dependence on 1 function 
of 1 variable. Thus for general class a; = 1 compatible PDE system 
we need three solvable Lie group to integrate it in quadratures: one 
group Gi of dimension equal to corang of the characteristic space n 
to perform the reduction {S,C£) — )■ (M, A), the second group G2 of 
dimension equal to corang of Aqo in M, and finally the third group G3 
of dimension rank(Aoo) — rank(A) (all actions should be transversal). 

Remark 2. A more general result is this: Consider a solvable Lie 
algebra G acting as transversal symmetries of A in M. Denote by 
A in M the quotient distribution. Then integral curves of A can be 
found from integral curves of A by quadratures. The number of involved 
integrals in the formula for the general integral curve is equal to the 
length of the derived series ofG. 

Let us consider an example from j Str] of a Monge equation S on 
y = y{x), z = z{x) with 3-dimensional solvable symmetry group: 

z' = z' + i:{z) + iy" + yf- (1) 

The group G = Sym(£^) is generated by the (prolongations of) vector 
fields cos X(9y, sin xSj^ on J°(M, M^) = M.^{x,y, z). 

If one (naively) substitutes y = h{x), then z{x) satisfies a Riccati 
equation, and so its solution cannot be found by quadratures (a similar 
problem occurs for general u = 1 class PDEs, so general reduction to 
ODEs from Theorem Rl does not necessarily yields a solution). 

The correct approach of Theorem [1] is to consider the quotient, i.e. 
to pass to the space of G- invariants M.'^{z, y" + y). A curve in this space 
is given by an equation y" + y = f{z). Substituting this back into ([1]) 
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we find the autonomous first order equation 

z' = z' + ^{z) + f{zf, 
which is easily integrable in quadratures. 

Remark 3. The previous naive argument uncover as follows. The 
curve in the plane Li^ = M.'^{z,y" + y) is specified via a parameter x: 
y" + y = h{x) and z{x) is given by z' = z"^ + ilj{z) + h{x)'^ . Since the last 
equation cannot be integrated in quadratures, the initial data (curve in 
LF' , a point over it determined by 3 constants) is not given explicitly, 
and so the result ceases to be given via an explicit formula. 

2. Integrable extensions and generalized symmetries 

Integrable extensions or coverings |KVj are mappings of PDEs £ — )■ 
B such that solutions of £ are obtained from those of £ by solving 
ODEs. For (underdetermined) ODEs the covering is just a submersion 
TT : (Af, A) — )■ (M, A), i.e. dxT^ ■ — )• A7r(x) is an isomorphism for any 
X e M. 

These coverings of systems of ODEs (or distributions) were studied 
in [AK] as they are useful in solving the system. Indeed a sequence of 
integrable extensions can decompose a given system into a sequence of 
1st order scalar ODEs. 

It is easy to see that quotient by the Cauchy characteristic of class 
u = 1 systems, which is basic for Theorem [Xj commutes with integrable 
extension. Thus it is enough to study integrable extensions of rank 2 
distributions. We will relate them to the symmetry approach of the 
previous section. 

For instance, we can write the symmetry reduction of Theorem [1] and 
Remark [2] via integrable extensions. Let p : M — )■ ]R(x) be function 
(projection) giving the independence condition. Write the equation for 
integral curves of A as F[x,u] = 0, where the latter is an ordinary (non- 
linear) differential operator and both F and u are mult i- dimensional. 

Then provided G has derived series G = Gi D Gi^i D ■ ■ ■ D Gq = 
with Abelian quotients Gi/Gi^i = ^ we can choose a coordinate Vi on 
Vi and have the equation for integral curves of A in this form: 

F[x,u]=0, v[ = Hi[x,u], = H2[x,u,vi], 

v'i_i = Hi^i[x,U,Vi, . . . ,f;_2], v[ = Hi[x,u,vi, . . .,Vi-i]. 

2.1. Existence of integral de-prolongations for (2,5) distribu- 
tions. Due to existence of normal forms (2, n) distributions have the 
structure of integrable extensions for n < 5. This holds true also in the 
first non-trivial case n = 5, where such distributions have moduli. 
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Theorem 2. A regular 2- distribution A on a manifold admits local 
submersion onto a 2- distribution in a ^-dimensional manifold (M'^, A). 

The claim foUows from (is equivalent to) a result due to Goursat: 

Theorem 3 ([Gsj, §76). A regular rank 2 distribution in a 5-dimensional 
manifold can be locally represented as the canonical distribution of the 
Monge equation £ : v' = f{x, u, u', u", v) . 

Indeed, the distribution of this equation £ C J^'^(]R, M^) is 

A = (9^. + Mi(9„ + M2<9„i + /<9^,<9„2), (2) 

which has the structure of integrable extension over J2(M,M), equipped 
with the canonical Cartan distribution {d^ + Uidu + <9mi , ) ! the 
projection tt : £^ — )■ is (x, m, Mi, U2, f ) ^ {x,u,ui,U2)- 

For completeness we give an alternative proof of Goursat's theorem 
(using vector fields approach instead of EDS methods). 

Proof Let the flag of the distributioE| be Ai = A, A2 = [A,Ai], 
A3 = [A, A2] = TM (we consider the general situation; in the other 
cases the distributions have normal forms and the statement follows). 

Consider the maps T : r(A) x r(A) ^ T{A*TM) and 9^ : r(A) x 
r(A) ^ r(A^TM) given by 

T(C,r^) = CAr^A[C,r7]A[C,[C,r/]], 

QoiCv) = rA[v, [C,v]], QiiCv) = T A [c, [C, [C,v]]i 

e^iCv) = T A [C, h [C,V]]1 03(C,^) = T A [r^, [v, [C,v]]]- 
A change of frame ( = a( + brj, fj = c( + drj induces the changes: 

r^ei(c, v) = a' ei(c, v) + a% (2G2(c, + ©3(^7, 0) 

+ afe^ (262(77, C) + ©3(C, V)) + 0i(^, C) + ^0o(r7, 0- 

where a = a ■ {arj + b(){b) — b ■ [arj + bQ{a). This implies existence of 
a solution f G C°°{M,RP'^) to Qi{C,v) = 0- 

Let us straighten ( = in a local chart M, and denote 

the quotient by = R^/C (i-e. M2 = const). Then the distribution 
becomes a U2-dependent vector field rj = A/C^in R'^. The Lie derivative 
L|j corresponds to the derivative by U2, which we denote by the prime. 



a b 
c d 



''By the commutator of two distributions we mean the distribution generated by 
the commutator of sections, e.g. r(A2) = [r(A),r(Ai)] etc. 
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Condition 6i(C,?7) = reads r] Arj' Arj" Ar]'" = 0, and we can assume 
the highest derivative can be resolved: 



By reparametrization of time U2 and scahng of rj we can achieve = 
ai = (in contrast the Laguerre- Forsyth canonical form). Then the 
equation is rj'" = a^rj" , and the solution \s rj = + U2ii + fd^, where 
fu2U2 7^ 0; ^''W'^v and ^o!^i ai'e M2-independent fields on = "K^/dy. 

Now in our general case the distribution (^o;^i) in is contact, so 
in the Darboux coordinates = dx + uidu, = du^. Thus we obtain 
local coordinates on M such that A has form ([2]). □ 

2.2. Non-existence of integral de-prolongations for (2, n) dis- 
tributions with n > 5. Dimensional count: generic rank 2 distribu- 
tion in M" depends on 2{n — 2) — n = n — A functions of n variables 
(quotient of sections of (2, ?T,)-Grassmanian by the pseudogroup of lo- 
cal diffeomorphisms) , while integrable extension depends on 1 function 
of n variables (in both cases: and some number of functions of fewer 
variables). Thus for n > 5 there are obstructions to existence of the 
structure of integrable extension over a lower-dimensional manifold. 

These obstructions are important relative differential invariants of 
the distribution. For example, in dimension 6 there are 2 relative in- 
variants vanishing of which characterizes possibility to represent the 
distribution as the Monge equation v' = f{x,u,u',u",u"',v)). 

It is interesting to notice that for n = 5 we get seemingly deter- 
mined system (the same functional dimension). And indeed, the one- 
dimensional distribution V has the property of de-prolongation (the 
vertical distribution of the projection vr) iff [V, A] C V + A. This writes 
as 4 equations on 4 functions specifying V (these latter can be taken as 
the 1st integrals, but then the system has order 2; it is better to write 
V via a generating vector field + Fidx + F2du + F^^dp + F^dq and take 
the components Fi as the unknowns). This 4x4 system is not however 
determined: it has every covector characteristic (direct calculation or 
the following observation - in the normal form of the previous section 
the totality of integrable de-prolongations has functional moduli - the 
general solution depends on a function of 5 variables). 

Remark 4. It is also interesting to try de-prolongation by rank 2 fo- 
liation, which indeed exists as a generic (2,4) distribution has Engel 
normal form, and so integrally de-prolongs to the contact (2,3) dis- 
tribution. The conditions for existence of such rank 2 distribution V 
are: 



III 



a2r]" + aiT]' + aoT] 



V 



[V, V] C V, [V, A]cV + A. 



(3) 
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This is a system of 4 equations on 3 unknowns (1st integrals ofV), but 
it is not overdetermined: again all covectors are characteristic! 

2.3. Generalized symmetries. A symmetry of the distribution A is 
a vector field ^ such that C^{A) = A. More generally a space Q of 
vector fields is a Lie symmetry algebra iff 

[Q,G]cg, [g,A]cA. 

Remark 5. Let us recall that a collection of differential operators (Fi) 
form a symmetry algebra for the PDF system S = {Hj = 0) iff 

{F^, Fp} = mod(Fi), {F,, H^} = mod{Hj), 

where {,} is the Jacobi bracket (we write the condition for simplicity 
in the case of scalar or square matrix equations), see [KLV] . 

// we are interested in compatibility of the systems (Ft) and {Hj), 
then the last condition changes to more general 

{F^,H^} = Omod{F,,Hj), 

see |KL] . Such F are called generalized symmetries, conditional sym- 
metries or auxiliary integrals. 

Basing on this remark we can treat distributions V satisfying condi- 
tion ([3]) as generalized symmetries. It allows the following symmetry 
reduction: if L C M is an integral curve of A, then in the union of 
l^-leaves meeting L the integral curves of A can be found via lower- 
dimensional determined ODE (for integral manifolds the corresponding 
criterion is a bit more complicated). If, in addition, V is obtained as the 
tangent distribution of solvable Lie algebra sheaf, the solutions can be 
found in quadratures. The generalized symmetries are more common 
than the classical ones. 

Example. Consider the symmetries of the Engel distribution, which 
is the Cartan (higher contact) distribution on J^(M, M) = M''(x, y, yi, y2)- 
In canonical coordinates it is A^ = (^i = + yidy + y2dyj^,C,2 = dy^)- 

By Lie-Backlund theorem the symmetries are lifts of contact fields 
on J^(]R, M), and so are defined by 1 function of 3 arguments. 

The generalized symmetries rj = dy + Xidx + X2dy^+^3dy^ (here unlike 
for symmetries we can normalize one of the coefficients by scaling) are 
defined by [^i f]] = mod A + {t]) which is equivalent to 

2/iAi - 1 l-yiM 1-2/iAi 

So the generalized symmetries depend on 1 function of 4 arguments Ai. 

It is often the case that a system (distribution) has no symmetries, 
but it admits generalized symmetries that can (partially) integrate A. 
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3. Integration of class u = 1 systems 

In this section we split the totahty of u = 1 systems into classes, 
and discuss transformations between them as a method of integration, 
r = dim H*'^{S) will be the total amount of PDEs in the systemj^ 

3.1. Type and complexity. We introduce the following rule for a 
choice of generators of the system S of class u = 1. Consider the 
orders of the PDEs in the system: fcmin = ki < ■ ■ ■ < kr = k^^x, which 
are taken with multiplicities irii = {j^j : kj = i} = dimH^~^'^{S). 

So the system £ is given by m^^ equations Fi i, . . . , -Fi of order ki, 
. . . , ruk, equations . . . , F^,™, of order kg = fcmax (s = r - nik, + 1). 

r 

We write S symbolically as ^ E^^ = Y^niiEi, and call the latter the 

i=l 

type of S. See [K] for the table of class a; = 1 systems of order fcmax < 5 
(this table works equally well for general non-linear systems). 

Let gi denote the symbols of £. Starting from some jet-level t the 
dimensions of these subspaces stabilize: dim^fj = 1 for i > t. This is 
equivalent to involutivity of the prolongation £;(*-'=max)_ 

Definition 1. Complexity of £ is the number x = Yli^oi^^^di ~ 

This number measures the amount of Cauchy data needed to specify 
a solution. It gives a partial order on the totality of class u = 1 systems. 
All our reductions will decrease the order. 

By definition all systems of class a; = (for ODEs the relevant 
complexity is the dimension of the solutions space) are taken to be of 
lower complexity than the systems of class u = 1. 

Lemma 2. Denote by £ the equation prolonged to the jet-level t = 
min{i : dim^fj = 1}; where it is involutive. Then dim£^ = x + t + 3. 

Proof. Since the base is 2-dimensional, we get dim£ = 2 + ^*^^ dimgi, 
whence the claim. □ 

Denote the Cartan distribution of £^ by C^^. We will assume that the 
system is prolonged to the level it is involutive, then dimC^ = 3. 

According to Lemma [1] the distribution Cs vcontains a Cauchy char- 
acteristic field ^ (the 1-dimensional distribution generated by it is 
unique). Denote by (M, A) be the (local) quotient by ^. This A is a 
rank 2 distribution describing the internal geometry of £. By Lemma 
|2]the manifold M has dimension /i = x + t + 2. 



'^Starting from this Section we restrict to base dimension n — 2. Familiarity with 
the Spencer cohomology H^'^{£) [S] is not cruciaL 
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3.2. Derived flags of a rank 2 distribution. Consider the strong 
derived flag of A deflned by Vi = A, [r(V,+i) = [r(Vi), r(Vi)] (where 
r(A) denotes the module of sections of the distribution A). The strong 
growth vector is the flnite sequence of dimensions (dimVj)[=;^, where 
r is the stabilization level (in the regularity assumptions, we adopt, all 
the ranks are constant). 

The weak derived flag is given by Ai = A, r(Ai+i) = [r(Aj), r(A)]. 
Notations 9^-^A = A,- are also used. The following cases are possible. 

I. The growth vector is (2, 3, 4, ... ). In this case by Cartan theorem 
C2I lAK] the system can be de-prolongec^, i.e. there exists another 



manifold M of dimension /i = /x — 1 equipped with rank 2 distribution 
A such that A = P(A) is the prolongation. 

The symmetries of A are preserved under passage to A, and the 
solutions are mapped forward in such a way that to any solution of A 
there corresponds a 1-dimensional family of integral curves of A. 

Thus passage to de-prolongation is a nice reduction of the system, 
for which the complexity x (it exists on both ODE and PDE levels) 
decreases. For linear class u = 1 systems this corresponds to the 
(generalized) Laplace transformation, see [K] and the next section. 

II. The distribution A is not completely non-holonomic, i.e. V,- 7^ 
TM. In this cas^ there are p = fi — rank Vr flrst integrals Ii, . . . ,Ip 
that pull-back to flrst integrals of the system £. We can flx the values 
of Ij and reduce the complexity of the system. 

For linear systems existence of such integrals means that the se- 
quence of Laplace transformations does not reduce £ to Ei but stops 
on a flnite type (class a; = 0) system [Kj. For non-linear systems the 
relative invariants that control existence of intermediate integrals can 
be calculated as generalized Laplace invariants of the linearization. 

III. The general case: the distribution is totally non-holonomic and 
not de-prolongable. Thus the growth vector is (2, 3, 5, . . . , /i). 

To flnd integral curves of A one can evolve Theorem [H or use inte- 
grable extension idea of Section[2]to decrease the complexity. Of course, 
due to calculations in Section [2l2] a generic rank 2 distribution on high- 
dimensional manifold has no integrable extensions, but distributions 
with symmetries do have such extensions. Indeed the symmetry reduc- 
tion gives the projection (with flbers almost everywhere transversal to 



^In [AKj the growth vector of the weak derived flag was considered. However 
this makes no difference at the first 3 elements of the sequence (2, 3, x, . . . ), where 
X = dim A3 can be 3, 4, 5. 

9A gain here it makes no difference if we consider weak or strong derived flag, 
only the length t can change. 
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A), that's why we can treat integrable extensions as generahzed sym- 
metries. Thus search of integrable extension is an integration method. 

Notice also that the quotient by Cauchy characteristic is not inter- 
changable with de-prolongations and restrictions to the level of first 
integrals, and it must be performed first. But sometimes the system 
needs to be prolonged for this. 

Example: Consider a compatible system £ C J'^(M) of type 2E^: 
Uxxx = F,Uxyy = G, with F,G e C°°{J'^M)). The system is not 
involutive on the level of 3rd jet^, so if we do not prolong the system, 
then the weak derived flag of the Cartan distribution Cs is (4, 7, 9, 10). 
In addition Cs has no Cauchy characteristics, while its derived d Cs has 
3 Cauchy characteristics, so that the pattern is wrong. 

The prolonged system £^^'^ C J^(M) is involutive and the reduction 
of Theorem |A] works - there is 1 Cauchy characteristic for the original 
distribution and one more for the derived. 

Consider for instance the system with F = jUy^y, G = \u'^yy 
The weak derived flag of the reduced (by Cauchy characteristic) sys- 
tem has growth (2,3,4,5,6,7,8,9), while the strong growth vector is 
(2, 3, 4, 6, 9). Thus there is one intermediate integral u^xy — \'^'^yy = c, 
and after de-prolongation both growth vectors are (2, 3, 5, 8) - the cor- 
responding graded nilpotent Lie algebra [H lAK] is free truncated. 

This shows importance of prolongation of the system S to the jet- 
level with dim^ffc = 1. 

3.3. Internal geometry of linear systems. Linear compatible PDE 
systems of class u = 1 were studied in [K] , which we briefly summarize. 

It was shown in that paper that such systems S (with dependent 
variable u) can be integrated via generalized Laplace transformation, 
which is a first order differential operator L : u v = Xu, with X 
having the same symbol as the characteristic vector field. 

Denote the system we obtain on the variable v hj S. It is also linear 
and compatible. Denote the inverse operator by L^^ : v ^ u. As 
proven in [K] only three different situations are possible: 

(1) S has class u = 1 and is a differential operator. 

(2) S has class u = 1 but L^^ is given by a finite type system. 

(3) S has class u = and L^^ is an integral operator. 

Case (1) is generic. If the itinerary of the transformations for £ meets 
only such equations, then Laplace transformation provides complete 
integration of the PDE system S. 



^^It has non-zero Spencer cohomology H^'^{£) = M, and the symbol is not stable: 
dim 53 = 2, dimgs+j = 1. 
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Moreover, under generalized Laplace transformation the complexity 
strictly decreases. Generically it decreases only by 1: x i— )■ x — 1. 

These results were obtained using the external geometry of £. Let 
us reformulate them in the internal language. 

Proposition 4. A generalized Laplace transformation for uj = 1 li- 
near systems is composed from the following maps in subsequent stages: 
some number of prolongations, a diffeomorphism, some number of de- 
prolongations. For an involutive system only two last steps are required. 

Proof. Indeed, from internal viewpoint the rank 2 distribution A, ob- 
tained from S via reduction by Cauchy characteristic, is a Goursat 
distribution (or Goursat in the leaves of the 1st integrals if the distri- 
bution is not totally non-holonomic). Since the Goursat distribution 
has the canonical normal form (see |Kul IMoj . we neglect singularities) 
the claim follows. □ 

Let us show how this works. We start with generic linear SE^ of class 
u = 1. Then in three Laplace transformations it becomes equation of 
type El (we refer to [K] for particular examples). We indicate the 
growth vector consisting of ascending by 1 integers, and indicate the 
internal coordinates on the equation: p,q,r,s,t are the classical nota- 
tions for the 1st and 2nd derivatives, and g is one of the 3rd derivatives). 



£ = 3E. 



S = Eo + E-i 



S = 2Eo 



£ = Ei 



\x,y,u,p,q 




x,y,u,p 




x,y,u 




x,y\ 


r,s,t,g 


>■ 


q,r,t,Q 


>■ 




>■ 


u,q\ 



TSdA 
(3,4,..., 9) 



(3,4,..., 8) 



(3,4,..., 9) 



(3,4,..., 6) 

t 

t 

rpg{2) 3 ^(2) 

(3,4,..., 8) 



T£d A 

(3,4) 
t 

t 

rpg{2) 3 ^(2) 

(3,4,..., 6) 



The first transformation is a diffeomorphism followed by a de-prolon- 
gation, the next one is a diffeomorphism followed by two de-prolonga- 
tions and the third one is of the same kind. 

If we choose non-generic SE^, then the route could be 3-^3 2E2 
El, so that the first stage contains more de-prolongations. 

Starting with 2E3 one needs to prolong once to follow the scheme. 
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Remark 6. Now we can explain decrease of complexity x via internal 
geometry. Since x = dim£^ — t — 3, where t is the order of involutivity, 
we see that x is defined correctly even if we prolong above the involu- 
tivity level (increase dim£^ and t equally). The difjeomorphism in the 
above proposition does not change the dimension, but it increases the 
order. Whence the claim. 

3.4. Closed form of the general solution. The shift along Cauchy 
characteristic is a characteristic symmetry for S, so its flow is tangent 
to solutions of this system (induces the trivial vector fleld on Sol(£^)). 
Thus S has closed form of the general solution iff the same is true for 
the reduced underdetermined ODE (encoded by A). 

Recall that closed form means possibility to represent general solu- 
tion parametrically through arbitrary functions and free parameters, 
which for PDE systems of class u = 1 writes 

(x,i/,«) = vl>(^,/(r),f(r),...,/(^)(r)). (4) 

Here is a M^-valued function with Jacobian of rank 3. For underde- 
termined ODE we shall remove one independent variable to the left, 
and the free variable <^ to the right. 

For rank 2 distributions the criterion for closed form description 
of the integral curves is known since E.Cartan. Namely in IC2I IKu] 
this is shown to be equivalent to A being Goursat, i.e. the canonical 
distribution C on the jet-space J'^(M, M), d = n — 2 = x + t. 

On the other hand we have demonstrated in [Kj that systems that are 
internally equivalent to ( J'^(M, M) x M, C x M) are internally linearizable 
and have no intermediate integrals. 

Intermediate integrals correspond to constants in the form of the 
general solution 

(x, y, u) = M/(^, /(r), /'(r), . . . , /(^)(r), Ci, . . . , c^). (5) 

The distribution can be transformed to Goursat-Frobenius normal form, 
namely it is internally equivalent to ( J'^-'"(]R, M) x M"' x M, C x x M), 
where m = codim(Aoo) is the codimension of the bracket-closure of A. 

Linearizability is not hampered by the additional constants. Thus 
we get the following statement. 

Theorem 5. General solution of a generic compatible system £ of class 
u = 1 and complexity x can be expressed in a closed form via a function 
f, its q < >c derivatives and x — q constants iff £ is linearizable by an 
internal transformation. 

Proof. Possibility to express solutions of a linear compatible systems 
of class w = 1 in closed form is proved in [K], so we need only to 
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demonstrate that q + m = x, where m is the amount of the first 
integrals (constants). 

The amount of derivatives to express all internal coordinates u„, 
|cr| < t, on £^ is — m. However the derivatives u„ are obtained from 
{x,y,u) via \a\ differentiations, so u shall be expressed in d — m — t 
derivatives of / only, and this number equals to x — m. □ 

3.5. Transformations of non-linear systems. Let us discuss some 
features of the transformations theory in the non-linear case. 

I. Quasi-linear systems allow some de-prolongations, but generic 
pure order k systems {kEf^) have none — after quotient by Cauchy 
characteristics the growth vector is (2, 3, 5, ... ). However if S is invo- 
lutive (t = fcmax) with different orders (/cmin < ^max), then we claim: 
The top equations are quasi-linear, and this implies existence of at least 
one de-prolongation. 

Indeed, provided that the characteristic is — Xdy, where A is a 
function on the jets of order /cmin, the top derivatives on the level k = 
fcmax must satisfy Ui^k-i = A*Mo,fc (this is due to the fact that ^ is the 
characteristic for all PDEs of £). Thus the PDE of order /Cmax in ^ can 
be chosen linear in top-derivatives. 

Consequently — XDy + pdu^ j. is the Cauchy characteristic of for 
some function p on fc-jets, and the two other generators of Cs are Vy 
and j. . A straightforward calculation yields that the latter field is a 
Cauchy characteristic for the derived distribution dCs, so the system 
can be de-prolonged. 

II. Re-covering the first integrals is the same as for the linear sys- 
tems. These restrictions introduce constants to the form of the general 
solution of similar as de-prolongations add derivatives to the form 
of the general solution. The two latter operations commute with each 
other and also with the projection of an integrable extension. 

III. If £ allows the structure of (i-dimensional integrable extension 
£ ^ £, then its solutions can be expressed via those of as m = L^{u), 
where is the resolution operator of a scalar ODE of order d. 

Provided a solvable Lie group of dimension d acts by transversal 
symmetries, the operator Lj can be expressed via the d-multiple quad- 
rature (I^T^ = J Ddr being quadrature by the parameter r). 

For instance, if A is maximally symmetric non-Goursat distribution, 
then its de-prolongation is fiat in the sense of Tanaka [TJ, and so has the 
structure of successive integral extensions over the rank 2 distribution 
in with G2 symmetry [AKj . So we get 
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Theorem 6. If £ has reduction A, which de-prolongs to a Tanaka-flat 
rank 2 distribution, then E can he solved in closed form and quadratures. 

Indeed, de-prolongations can be interpreted as non-linear Laplace 
transformations with differential inverses (this yields a closed form over 
the solutions of the reduced system A), while Tanaka flat rank 2 dis- 
tributions A project via integrable extensions to the Hilbert-Cartan 
equation |AK] (so its integral curves are given in quadratures). 

Thus we get the next easy case (after linearizable systems) of explic- 
itly integrable class a; = 1 systems, which are reduced to symmetric 
Monge systems (these latter were classified in |AK] ) . 

Remark 7. // S is an integrable extension, then its solution form can 
be specified via particular quadratures, like it is done in \Ci\ IG3I in the 
case of rank 2 distribution in dimension 5. 

Finally let us mention that the solutions obtained via the proposed 
method are usually not of Moutard type, where the latter means that 
u is expressed directly as a function of x, y (in the above approach all 
X, y, u are expressed as functions of the additional parameters r). 

Indeed in order for the discussed transformations to be Moutard, 
they shall preserve the vertical fibrations, i.e. the Cauchy character- 
istics for the derived systems and the generalized symmetries must be 
tangent to the vertical fibers (the 1st integrals are always independent 
of the base coordinates x,y). 



4. Examples of symmetric PDEs 

4.1. Model reductions to ODEs. Consider the following compatible 
class u = 1 systems of the type kEk'. 



AEi : u, 



'XXXX 



3^3 : u, 



A, u.. 



2E2 : u.. 



'XXX 



•xxxy 



A, M; 




'XX 



•xxy 



A, U: 



'Xxyy 




yyyy 



A^ 
4 



= — etc. 
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The reduced growth vector£l| and the generators of the weak derived 
flags are the following: 

(2,1,2) : (61,61,62,63,63), 

(2, 1, 2, 3) : (61, 61, 62, 63, 63, 64, 64, 64), 

(2,1,2,3,4): (61,61,62,63,63,64,64,64,65,65,65,65") etc. 

Commutators are given by [ei,6'J = 62, [61,62] = 63, [e'1,62] = 63, 
[61,63] = 64, [61,63] = [61,63] = 64, [61,63] = 64 etc (the commutators 
of 62 and 63 and others are zero) — these yields the structure of graded 
nilpotent Carnot algebra associated to the weak derived flag [Tl lAKj . 
The corresponding Monge underdetermined systems of ODEs are: 

(2,1,2): y' = \{z")\ 
(2,1,2,3): y" = \{z"')\ u' = \{z"')\ 

(2,1,2,3,4): /' = l(^-)2, n"=|(0^ t;' = l(/Y etc. 

This follows from the explicit form of the generators. Indeed, let us 
demonstrate this, for simplicity, in the case 3i?3. 



61 = -V^ = -(5^ + u^du + u^xdu. + Uxyduy + \du,^ + ^du^y + ^duyy), 
e'l = dt, 62 = (9„_ + Xdu^y + X'^duyy, 

63 = + Xduy, 63 = du,y + 2Xduyy, 

64 ~ du, 64 = c^My, 64 = Idr^yy. 

In this list we have omitted the generator of (rank 3) Cartan distri- 
bution T>y because the Cauchy characteristic equals T>y — XD^ and we 
need to quotient by it. 

Now to perform the quotient one has either to pass to invariants, 
or to restrict to the transversal of the action. We choose the second 
approach: Change the notations z := u, z' := u^, z" := u^x-, z'" := A = 
Uxxxi y '■= Uy, y' := Uxy, u := Uyy and we are done. 

Remark 8. To recover the solutions of S from the solutions of the 
ODE system, the first method must be used. For example, the first of 



the equations — system 2E2 from Ci — has the following invariants 



of shifts along Cauchy characteristic (uxx = X) 

t = —A, w = —2u + 2yUy + 2xUx — x'^X — xyX^ — |y^A^, 
z = Uy — UxX + |xA^ + \yX^ , zi = Ux — xX — \yX^, Z2 = x + yX. 



^-"^These are the ones we have used in |AK| : we pass from the usual growth vector 

(ni, 712,71-3, . . . ) to (ni,n2 - Til, 713 - 7l2 . . . ). 
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Thus we can express the general solution parametrically as 

X = z"{t) + st, y = s, u = SZ + z'z" - - ^t z"^ - \t^s z" - \t^s^, 

where the two functions z = z{t),w = w{t) are related by the Hilbert- 
Cartan equation w' = {z")"^. 



Another interesting sequence of equations, considered in |Gi| , is pro- 
vided by 

2n 



u 



xy 

X + y 



V^xUy. (6) 



These PDEs are Darboux integrable with intermediate integrals of or- 
der (n + 1), see |XF] . 

In order to integrate ([6]) let us linearize it via Goursat substitution 



(j2| p = a/m^, q = y/u^, which leads to the system 



n n 
Py = — : — Qx = — —P- (7) 

X + y X + y 

Then P^-intermediate integral can be found as an ODE on q of the 
form (I?;,,-intermediate integral is obtained similarly via an ODE on p) 

n 

L = y —^l^q, = 0, (8) 

where = T>yq. The condition of intermediate integral — = 
on (171) — is an overdetermined linear system on constants ctj. With 
normalization = 1 its unique solution is given by 

n^{n - If ■ ■ ■ {n - i + If 



a 



n—i -1 



Consider now the overdetermined compatible system ([6]) + ([8]). It has 
class a; = 1 and type E2 + E^, m = n + 1. The Cauchy characteristic 
is Vy + (pmdg^ for a properly chosen function yj^. 

Reduction along Cauchy characteristic (which again can be inter- 
preted as intersection with the level of the Pa.-intermediate integral) 
yields a rank 2 distribution on a manifold S of dimension (2n + 3). This 
system has n de-prolongations, and so it reduces to a rank 2 distribu- 
tion on {n + 3)-dimensional manifold. 

The symmetry analysis (done by I.Anderson) coupled together with 
unique symmetry model for rank 2 distributions |DZt lAK] implies that 
this rank 2 distribution corresponds to the Cartan distribution of the 
Monge equation 

y' = {z^""^?- 
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In particular, we recover the result (known to Goursat) that the so- 
lutions of equation ([S])n=2 are expressed via the solutions of Hilbert- 
Cartan equation. 

Remark 9. Thus we realize two boundary lines from the Zoo of types 
in [K] — the bottom and the diagonal — as the most symmetric PDEs 
in its class (both types of PDF and the reduction are fixed). 

Here we refer to the contact symmetries (internal symmetry group 
of class u = 1 systems is infinite- dimensional), which turns out to be 
isomorphic to the internal symmetry group of the reduction by Cauchy 
characteristic in most cases. This Lie-Bdcklund type theorem will be 
discussed in details elsewhere. 

4.2. Representation of solutions. Linear class a; = 1 compatible 
systems £ have their solutions expressible in Moutard form, i.e. u is 
a parametrized function of {x,y). In general, the closed form solution 
force all of the variables (x, y, u) to be expressed via parameters (C, t) 
— this is equivalent to linearizability, see §3.41 

Linearization map does not commute, in general, with the projection 
to the base {x,y), that's why many exactly solvable system do not 
possess Moutard form. We illustrate this with two 2nd order examples. 

Example 1. Consider the system ([6]) + (|8])„=i: 

^xx 2 ■ , Uxy 2 - . (9) 

X + y X + y 

The reduced growth vector is (3, 1, 1, 1), so S is internally linearizable 
and is solvable via generalized Laplace transformations [K]. It is how- 
ever non-Moutard. 

To see this let us describe the general solution. Goursat substitution 

Ux = Uy = (10) 

linearizes the equation 

p p q 

Px = — , Qx = — , Py = ■ • (11) 

X + y X + y X + y 

Notice that this has vector growth (3, 1, 1), and so one could suggest 
it serves as de-prolongation of Qj, but it does not. The reason (as 
shown below) is that it is impossible to de-prolong preserving the base 
coordinates x, y. In fact, (IHl) is an integrable extension of ( ITT]) via ( |T0|) . 

The next step is to observe that the last equation of ( ITT]) can be 
used as a definition of q and the second equation is then the differential 
corollary of the first. Thus we can restrict to the equation Ei. 

Px = (12) 
X + y 
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with the reduced growth vector (3,1). Inverse transformations are: 
integrable extension 

Ux = Uy = {x + yfpl 

to the space M^(a;, m, Ma;, My) and then prolongation to the original 
equation £ = M.^{x, y, u, Ux, Uy, Uyy). 

Let us return to the closed form of the solution. Solving (fT2l) we get 
p = ip{y)/{x + y). Integrating the Frobenius system 

[x + yY V x + yJ 

yields the solution 

u = <P{y) - ^ with <P'{y) = i^'iyf. (13) 
X + y 

This latter constraint is internally equivalent to the Engel distribution 
(or to J^(R, R)) and the equivalence is given explicitly by 

y = a"(r), i> = Ta"{T)-a'{T), ^' = t, = rV"(r) - 2r(T'(r) +2(T(r). 

Thus alignment of the equation £ ^ to the jet-space J°'^(R, R^) = 
M^(x, r, a, (Ti, (72, da) is the following: 

X = X, y = o"2, u = T a2 — 2to\ + 2cr - 



X + (T2 



_ (ra2 - (y\f _ (rx + Oxf 
[x + [x + 

\2 



Uyy = 2(tx + ai] 



{{x + - {rx + cTi)a3) 



-v--^ (a: + a2)V3 



This determines a diffeomorphism J'^'^(R, R^) — )■ £. The inverse map 
£ — >■ J'''^(R, R^) is also given in differential-algebraic form 



X = X, T = y/Ux + O = -{U - XUx + yUy - 2Xy/u~lL^), 

2{x + y)^. 

(^1 = y^nhl - X^JUx, (72 = y, 0-3 



Thus Laplace transformation £ = 2E2 ^ Ei = {^^ = {f\ has the form 



2Uy + (X + y)Uyy ^ 

dx 



1 



X = X, r = ^Ju^ + ^/u^, Cr = 2 ~ •^''^^ ~^ y^y ~ 2Xy/UxU 



and it decomposes internally into the composition of the diffeomor- 
phism £ — )• J°'^(R, R^) followed by the double de-prolongation (we do 
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not change the original Cauchy characteristic, which is the direction of 
the first factor in J0'3(R,R2) = R x J3(R,R)): 

J°'^(R,R2) ^ J°'^(R,R^) =R^{x,T,a,ai). 

Let us explain why 2E2 (E]) cannot be transformed to Ei by a Laplace 
transformation preserving the {x,y)-ha.se (Moutard type). 

We wish to find a relation on x,y,u,Ux,Uy excluding the above func- 
tions (j){y),tlj{y). But this is impossible since 

yju^ = — ■ — , y/uy = ip [y) ■ — 

X + y ^ X + y 

and so u,Ux,Uy are algebraically independent. 

Another approach is to show that the constraint (f)'{y) = ip'ivY i'^ ^^e 
form f[T^ is equivalent to the standard Engel distribution on J^(R, R) 
internally, but not externally. Indeed, no point transformation can 
map the above constraint to the equation = 0, since their point 
symmetry groups have dimensions 10 and 00 respectively. 

Example 2. Another interesting system, discussed in §4.H is the 
Cartan involutive 2E2 model 

^xx -^5 ^xy ' '^yy 3^ ' (■^^) 

The reduced growth vector of its (5-dimensional) reduction £ is (2, 1, 2), 
so Laplace transformation in the sense of linear theory does not exist. 
But £ has the structure of integrable extension over Ei, namely over 
the gas dynamics equation 

Vy = vvx. (15) 

Indeed, this latter is just the compatibility condition on the parameter 
\ = V along the Cauchy characteristic. The transformation from f[T^ 
to f[T^ is a composition of a 3-dimensional integrable extension and 
the de-prolongation: 

n prol y / ext . 

R {X,y,U,Ur„Uy,Uyy) " " "> R {x , y , U , , Uy , Uyy , Uyyy) " " R {x,y,V,Vy). 

Since (fT5|) is clearly not of Moutard type, this explains that ( IT^ is not 
of Moutard type. Its solutions are though expressible via the solutions 
of Hilbert-Cartan equation. 

This makes fll4p internally equivalent to de-prolongation of the equa- 
tion ([S]) + (IHDn=2 of type E2 + E3 considered in §4.11 Another equation 
equivalent to f|T^ is thus the compatible 2E2: 

W^X = 0, + XW^y -Wy = Q. 
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4.3. Other symmetric models. In [Cij Cartan considers also sub- 
maximal symmetric systems £: 



2 



Tfl 

2E, : u,, = A, u,y = A™, uyy = ^A^"^-\ (16) 

Its contact symmetry algebra has dimension 7, and this is the next 
possible number after the maximal finite value 14 for dimSym(£^). 

The Cauchy characteristic is ^ = Vy — m\^~^Vx, and the reduced 
system £ can be found by restricting to the transversal y = const to ^. 
In other words, the rank 2 distribution of £ is given by its generators 

which after a change of coordinates is identical with the Cartan dis- 
tribution of the Monge equation (m ^ 2, 1 - the exceptional cases 
corresponding to dim Sym(A) = 14) 

w' = {v'T- (17) 
The higher analogs of (1161) are straightforward: 

2 3 

■ u =a u =a'^ u = a^™"^ u = a^'"-^ 

. Lixxx "'xxy " ; xyy 2,771 1 ' ^^"^ SwT, 2 

Dimension of the contact symmetry algebra here is 10 (this is readily 
checked with the help of Differential Geometry package of Maple), 
while dimension of the maximal symmetric non-linear model (system 
3^3 from SI]) is 12. 

The reduction by Cauchy characteristic is the following underdeter- 
mined ODE system: 

^ ' ' 2m- V ' 

It is a (3-dimensional) integrable extension of Monge equation (fT7|) . 

Similarly we construct equations AE^ etc. The sub-maximal sym- 
metric class CO = 1 system of type nEn is 

c Tn'^ 

-<V+^ : l<i<n 



I ' im — i + 1 

Its contact symmetry algebra has dimension |n(n + 1) + 4 (vs. the 
maximal dimension ln(n + 1) + 6 — the details on this calculation will 
be presented elsewhere). 

It is an interesting open problem what are the sub-maximal sym- 
metric PDE systems of the type E2 + En and what are (sub-)maximal 
models for the other types from the Zoo of [K]. 
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